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A unique polar nephelometer was designed and constructed for the 
measurement of atmospheric particulate characteristics. The nephetometer 
produces visible light from a self-contained laser to irradiate an air 
sample drawn into the instrument. The light scattered from the particu- 
lates and molecules in the sample is detected as a function of scattering 
angle for each of four different incident 1 ight polarizations. These 
measurements are used to determine the particulate scattering matrix 
which is a function of the size, shape, and index of refraction of the 
particles. The region of sensitivity for the measurements corresponds 
to the size range of particles that strongly affects visible radiative 
transfer in the atmosphere, which is the primary application for the 


derived information. 
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I. Introduction 

Atmospheric aerosol size distributions have been shown to be signifi- 
cant in numerous fields of study. Atmospheric particulates are essential 
to cloud formation, and clouds have a dominant effect on radiative transfer 
in the atmosphere. 1 Even for clear skies, the particulates affect the 
heating or cooling of the atmosphere through scatter and absorption of 
energy. Subsequently, the solar energy received at the surface and infra- 
red emission into space are also affected by aerosols. These combined 
factors enter into the analysis of the consequences of particulates in the 
atmosphere. Since only minimal data are available on atmospheric aerosols, 
the prediction of aerosol impact on climate is severely limited. 

Due to the importance of aerosols and the broad range of sizes 
occurring naturally and artificially, numerous techniques have been 
developed to measure particulate characteristics. Direct sampling is the 
most straightforward technique, but low counts, biased samples, and slow 
analysis leave much to be desired. Nuclepore filters 2 have been used 
for sizing the smaller particles (-below 0.1 pm in radius). Knollenburg 
counters 3 1 ikewise measure physical size. The technique involves the 
occultation by particles of a collimated light source projected on an 
array of photodetectors. The Knollenburg counter has a dynamic range of 
about one order of magnitude that Is set for the size range of interest 
down to 1 pm. Electrostatic devices such as the Whitby counter 4 utilize 
particulate mobility as a means of sizing aerosols. The Whitby has a 
usable range down to about 0.005 pm in radius. The maximum size limit 
is about 0.1 pm for relatively clean air such as in Tucson. 
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The passage of light through a medium produces scattering and absorp- 
tion which ts dependent on the medium’s properties. Particulates scatter 
tight according to their size, shape, index of refraction, and the wave- 
length, intensity and polarization of the Incident ltght. The information 
encoded by the aerosol on the scattered light is a major source for aerosol 
analysis of particles larger than 0.1 pm. Multiwavelength radiometers 
which measure the extinction of the direct solar beam and solar aureole 
detectors which measure near forward scattered light have been used to 
obtain aerosol size distributions which are averaged over the entire 
depth of the atmosphere. 5 

Within the optically active region (where atmospheric aerosols 
scatter light appreciably), the size distribution and refractive index 
are basic parameters which can be inferred with a polar nephelometer. 
Typically, the terminology "nephelometer" encompasses those devices 
which determine particle size and/or concentration with scattered light. 

They are small-scale devices that incorporate their own light source. 

One group includes those instruments which have a wide viewing angle set 
at a fixed position. Both particle counters 6 and integrating nephelometers 7 
fall into this category. A second group consists of polar nephelometers, 
which have a small field of view of scatter, but the angle at which the 
scatter is observed may be changed. These Instruments typically measure 
the scatter from a single particle or from a volume element that includes 
many particles. 8 * 14 

Bistatic lidar 15 is essentially the same principle as the volume-type 
nephelometer. The measurements are taken, however, from scattering volumes 
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which may ba kilometers away. Lldar has the desirable feature of being 
able to measure a<***oso1s at various heights In the atmosphere, but effects 
such as extinction due to the long light paths and background noise must 
be dealt with. Also, the ability to measure scattering In the forward 
direction is limited due to the geometry of the system. 

This paper's research encompasses the design and construction of a 
volume-type polar nepheiometer that will determine the foor elements of 
the MIe scattering matrix for particulates. 

II. Scattering Theory 

To evaluate a broad size distribution, it is necessary to sample a 
number of particles in an unbiased manner. The approach here is to 
measure the scattered light in situ from an air volume large enough to con- 
tain a statistically significant sample of particles over the range of 
sizes where scattering is appreciable. The scattering can be divided 
into two separate components for theoretical considerations, molecular 
and particulate. 

A. Molecular Scattering 

First consider the scattering from diatomic air molecules, which 
Is well-known from Rayleigh theory. The scattering equation can be 
written as 

I - FI, 

— s i 

where £ is the appropriate four-by-four scattering matrix and £ has the 
units of power per unit area per unit solid angle with the underlined 
notation indicating a matrix. I represents the scattered light vector, 
and I_j designates the incident light vector. 
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The general scattering matrix for the Rayleigh case as presented by 


van de Hulst 16 Is substituted Into Eq. (1) yielding 
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The incident and scattered light vectors are presented In their standard 
Stokes polarization notation. Ij and I* are the parallel and perpendicular 
radiances of the beam referenced to the plane that contains the incident 
and scattered light vectors. U Is a measure of the plane of polarization, 
and V describes the ellipticity of the beam. An alternate representation 
uses I ■ Ij + 1“ for the first element and Q ■ Ij - I* for the second, 
r ts the distance to the scatterer from the point of detection, X is the 
wavelength of the incident light, 0 is the scattering angle (see Fig. 1), 
and A and B are parameters which arise due to polarizab' 1 1 ty of the non- 
spherical molecules. The nonspherical effect Is dealt with by averaging 
over all possible orientations. Averaging applies quite well to the cases 
for nitrogen or air near standard temperature and pressure since very large 
numbers of molecules are In even the minimum scattering volume of the instru- 
mental arrangement. 

The depolarization ratio, p ■ may b e substituted into Eq. (2) 


giving 
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This Is easily transformed [with F' representing the four-by-four matrix 


in Eq. (3)] to 
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where m is the Index of refraction for air and N is the molecular concen- 
tration. Details of this derivation are found elsewhere. j^ e jQyy 
format is more convenient for this work; therefore replace F' with 
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where Q, indicates that the alternate terminology is in use. 1^ is changed 
from a per molecule basis in Eq. (4) to a per unit volume (V) notation, 
giving 
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where the t subscript Indicates standard values. Finally, pressure (p) 
and temperature (T) corrections are necessary to obtain 


2i»ir 2 (m -l) 2 VpT 

h* ' A T r r N t (S-Tp')p^T f-'d —I (7) 

B. Aerosol Scattering 

The aerosol contribution to the scattered light is complicated because 
the shape of an aerosol may be highly irregular and its composition may 
be inhomogeneous. These factors make it difficult to theoretically cal* 
culate the scattered light from even a single particle. By making the 
appropriate assumptions, such as that an effective refractive index holds 
and that the particles behave as spheres, the tractable Hie theory can be 
applied to predict scattering from a collection of particles. 13 It should 
be noted that for radiative transfer purposes it is not essential that the 
particles be modeled correctly. Hot/ever, it Is necessary that the model 
produce the same effect on the incident radiation field as the actual parti- 
cles do; thus, the reasoning behind such terminology as effective index of 
refraction. Continuing with the spherical case, one finds 
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where the matrix elements are defined by van de Hulst. 16 
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Expansion In terms of the amplitude functions, S, for Fq yields 
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Designate the bracketed portion of Eq. (9) as (where the subscript refers 
to the aerosol), end the scattered radiance for the aerosol becomes 
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which has Eq. (6) as its molecular counterpart. 

C. Total Scatter 

The complete theoretical representation of the scattered light is 
found by combining the molecular and aerosol parts giving 
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Due to the averaging and symmetry considerations made earlier, the 
theoretical scattering Is independent of the azimuthal angle, and Is 
only a function of the scattering angle (see Fig. 1). Even more generally, 
one can observe that the Mie scattering angular variation is similar to the 
Rayleigh pattern for particles whose radius is somewhat smaller than the 
incident light's wavelength. Unfortunately, this means that there is no 
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slz« Information from small particle angular scattei tng. Larger particles 
(Fig. 2) do have distinct scattering patterns, however, which make sizing 
possible as varies wfth both particle size and angle. 

III. Nephelometer Design 

Mie theory predicts some fraction of the incident energy at a 
specific wavelength for a given particle size and refractive index is 
scattered into a given direction. Due to previously discussed symmetry, 
a single 0 to 1 80 degree sweep totally characterizes a particulate's 
scattering pattern. Accordingly, if the scattered light is measured at a 
series of angles, information could be obtained about the particle's index 
of refraction and size. 

Previously, it was shown that there are four unknowns in the Mie 
scattering matrix. Since each measurement of the scattered radiance gives 
a single equation with four unknowns, four independent measurements allow 
the solution of the scattering matrix at a given angle. 

Constder the polar nephelometer which measures scattered light as a 
function of angle. Outside air with its particles is continuously drawn 
through the enclosed system by an exhaust fan (Fig. 3). In order to keep 
the nephelometer relatively compact, a folded light path Is used from the 
laser to the scattering region. The 0.51 ^*5 ym line of a pulsed argon laser 
provides the light source, so that essentially a single wavelength is 
scattered from the air. A photomultiplier tube (PMT) is used to detect 
the radiance of light scattered by the molecules and aerosols in the field 
of view of the receiver optics at a given scattering angle (Fig. A). This 
allows sufficient gain so that the weak scattering is above the noise level. 
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The scattering angle may be varied over a wide range (2 to 178 degrees), , 
but at any (me setting, only a narrow range of ±.25 degrees of scattering 
Is observed. The remaining variables are physical parameters of the aero- 
sol Itself. 

For accurate scattering measurements, a number of calibration factors 
must be considered. The PMT integrator Is gated on the laser pulse to 
avoid integrating noise during the off cycle of the laser. As many pulses 
as desired may be used to Improve the signal to noise ratio since the 
voltage on the integrator Is proportional to the light scettered onto 
the PMT. An LED acts as a stable light source for calibration of the PMT 
gain which varies (particularly for a change in supply voltage). The LED 
ts mounted in optical cm tact with the PMT and produces a pulse between 
laser pulses which Is gated and integrated in a fashion similar to the 
scattered tight signal. The LED Is temperature-controlled at 50*C with a 
heater, so that its output will remain virtually constant. The laser power 
monitor (PM) signal is likewise integrated to that any changes in energy 
Incident on the scattering v-’ume can be determined. (Attentuation by the 
scatterers is assumed negligible for purposes of thts measurement.) Since 
this monitor signal is large, a silicone photodiode, which does not experi- 
ence the same variability of sensitivity as a PMT, is used as a detector. 
Noise is effectively cancelled by zeroing the integrators with the laser 
light blocked out. 

Next, the measured scattered light is corrected for changes in the 
scattering volume with angle, in sensitivity of the detector with distance 
of angle off axis, and in laser energy axially. The correction curve for 
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that* combined factors closely follows sec 9. If P(0), the scattering 
phase function, varies significantly over a given scattering volume, the 
problem Is not as straightforward. Since the most rapid changes In P(6) 
normally occur near forward scatter, a check Is made by finding If 
£ VjPfej) 

freTTT:-' ' ,2) 

J J 

at those angles (where Vj Is the laser energy times the detector sens I* 

tlvity for the volume element , J). When Eq. (12) holds, then P(6) £ v. 

.1 J 

may be used as an adequate approximation where p(e) Is calculated at the 
angle of the central volume element. The specifics of the scattering 
geometry require careful attention and are covered in a rigorous manner in 
Hansen. 16 

The light leaves the laser linearly polarized. A polarization rotator 
is used to change the angle of linear polarization to any desired value. 

A 1 /4-wave plate may be inserted into the optical path to change the polari- 
zation to circular. Values of polarization rotator angle and quarter-wave 
plate angle are chosen to give four independent incident light polarizations 
and therefore allow a unique solution for the aerosol scattering matrix 
elements. Subsequently, some light scatters from a sample volume and 
passes through a polarizer which Is set at k$ degrees to the scattering 
plane. (The polarizer is used in front of the detector optics so that 
Djj and Sji appear In the radiance measurement.) The detector then 
gives a reading of the transmitted energy. The PMT integrator voltage 
(corrected for the discussed factor?.) ts proportional to 0 P(8), the 


volume scattering cross section times the phase function matrix. The 
values are different for aerosols and molecules so using a constant G 
which accounts for the terms In common gives 


13 


PMT - 



f 

1 

e 

polarizer 


matrix 


l 

J 




oj) 


The PMT voltage corresponds only to the radiance magnitude element, X, 
which Is indicated by the brackets* Pure nitrogen gas can be introduced 
Into the nephelometer. With only nitrogen for a scetterer,. the value of 
G may be solved as P^ and are determined from Rayleigh theory and 
Is aero. For the case of aerosols and molecules, the value of 6 ^P^ 
may be determined once G Is known. 

In summary, this instrument evaluates the aerosol characteristics 
through the determination of the scattering matrix at a given angle. The 
matrix is measured by automatically cycling the polarization of the laser 
light source through four different settings at the scattering angle. In 
the automated operation mode, measurements at successive scattering angles 
separated by I, 2, k, 8, or 16 degrees are made starting with 178 degrees 
scattering (Fig. 5). 

During the operation of the nephelometer, the integrator voltages and 
scattering angle are digitized by a four and one-half digit autoranging 
multichannel voltmeter, then punched into paper tape. The tape is used to 
create either a card deck or magnetic tape file. The data are then read 
into a computer program which solves for the four aerosol matrix elements 
for each angle by simultaneous equations with the four incident light 
forms. Some other distinct characteristics of the system follow: 
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1. Only one polarization is incident on the PHT to eliminate 
polarization sensitivity. 

2. A temperature-controlled LED is used as a gain monitor for the 
PHT. 

3. The system is automated for data acquisition. 

4. A small detector unit allows easy angular rotation for changing the 
scattering angle. 

5. An effective light trap (made of flat polished absorbing glass) 
for the unscattered part of the laser beam allows large scattering 
angle measurements. 

6. A pulsed laser is used so that peak power is high for a good S/N 
ratio while the laser itself is small and air-cooled. 

IV. Calibration 

To verify the proper operation of the system and to calibrate the 
nephelometer , the system was sealed and flushed with filtered air, then 
nitrogen gas. This procedure lowered particulate counts in the nephelometer 
to below the measurable level of an expans ion- type counter (less than 30 
particles per cubic centimeter) . Measurements of scattering that were sub- 
sequently taken showed reasonable agreement with the expected molecular 
scattering (Figs. 6 and 7)* Errors of 1 to 2 percent were observed at 
angles near forward and backscatter. These increased to 5 to 10 percent 
near 90 degrees. These values are felt to be sufficiently accurate as mul- 
tiple measurements can average out the error significantly for calibration 
purposes. Also, since the problem is one of signal to noise and since 
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aerosol scattering Is significantly greater In the forward direction which 
Is of primary Interest, Instrumental errors less than 2 percent would 
be obtained. The absolute magnitude of aerosol measurements Is determined 
by the system calibration constant obtained with the known magnitude of 
molecular scattering from theoretical considerations. 

V. Results 

Data are presented in Figs. 8-11 for particulate scattering matrix 
elements measured in Tucson, Arizona. The sampling location was about 
10 meters above ground level. Measurements were typically made In the 
evening when the atmospheric conditions were more stable in order to 
minimize changes with time in the particulate populations. Measurement 
time was approximately 1 minute for each scattering angle. An expansion- 
type particle counter and repetition of nephelometer data were used to 
verify aerosol stability. Observed particle counts were typically 
20,000 per cm 3 , and periods of stability lasted 15 to 30 minutes. 

On examination, the elements were generally observed to behave in 
a manner consistent with expectations. A limited amount of experimental 
data of a similar nature have been publ tshed. 11 * 13,2 ° Grams et al. com- 
pares most favorably as they also used atmospheric aerosols, but detailed 
quantitative comparisons are not valid unless simultaneous measurements 
are used. Curve fits are presented as a means of providing an analytical 
representation that follows the data reasonably well. 

The Mj and M2 elements (Figs. 8 and 9) are fitted with a double 
Henyey-Greenstein function using two scaling constants. The Henyey- 
Greenstein function is appropriate to use for these elements due to 
the phase function characteristics associated with both the data and 
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the Henyey-Greenstein function. The double function Is necessary to model 
both the sharp forward peak and the relatively flat backscatter region 
closely. An alternating gradient correction, then Gauss-Newton correction 
method Is used to converge to a best least squares fit. Error bars are 
determined by a percentage of the matrix element value based on the magni- 
tude of the radiance measurement which Is strongly angle dependent. A 
minimum absolute error Is also applied for very low values. 

For the S 2J and 0 2 j elements (Figs. 10 and 11), a log formula was 
used for fitting the data. Sign changes were handled by taking the log 
of absolute values of the data, then renewing the original sign. With 
the data in this form, a sixth order polynomial fit was calculated. Some 
of the S 2 j and D 2 j values are questionable because fast changes of sign 
or magnitude with scattering angle are not expected for a size distribution 
which normally has a smoothing effect. However, these anomalies occur for 
low values which have the greatest percentage error. Also, the statistics 
for large particles are the worst due to their low numbers and large cross 
sections. Subsequently, variation of the large particles in the scattering 
volume (particularly near 90 degrees scatter) might reasonably be the 
cause of some irregularities. 

VI. Error Analysis 

Errors in the determination of size distribution from nephelometer 
data fall into three major categories. Sampling errors are caused mainly 
by obtaining a biased sample of particles. This might happen through 
losses ?n the sampling tube which channels outside air to the nephelometer. 
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However, the sampling tube Is 8 centimeters In diameter to minimize 
losses. The effect due to changes In the number of very large particles 
in the scattering volume over the measurements In a set may also produce 
a bias. Forward scattering angles, where the large particles scatter 
strongly, have the larger scattering volumes which tend to offset this effect. 

A second category of errors is designated instrumentation errors. 
Scattering angle error is typically lfess than ±0.1 degree between the 
axial ray of the detector optics with respect to the axis of the laser 
beam. A detector acceptance angle of ±.25 degree causes an additional 
error due to a range of accepted light scattering angles. These factors 
become most critical very near forward and backscatter where scattering 
radiance changes with angle are greatest. Errors due to optical train 
misalignment and PHT nonlinearities are more difficult to assess. Measure- 
ments of the PMT and integrators show errors of less than 1 percent for all 
except minimal measurements such as occur with molecular scatter at 90 
degrees. Optics alignment Is somewhat difficult to maintain. Although 
adjustments are made to better than 1 percent, several percent is some- 
times measured at subsequent times. 

Finally, data analysis itself has inherent errors. Quadrature error 
and other factors due to data truncation are negligible compared to other 
errors. The main problem is the error due to approximating real aerosols 
as Mie particles. This effect could easily be 10 percent and dominate 
the other errors, but no evaluation is attempted here. 
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VII. Summary and Conclusion 

A polar nephelometer with a number of unique features has been developed 
for use tn studying atmospheric aerosols. The nephelometer has sufficient 
sensitivity to detect molecular scatter from air. This allows the instru- 
ment to measure scattering from very clean air and to use pure molecular 
scattering for calibration. A relatively compact system is achieved by using 
a folded light path with an air-cooled argon laser for the light source. 

A small, sensitive detector unit permits easy angular rotation for changing 
the scattering angle. A narrow detector field of view of +\/k degree of 
scattering Is used along with a single wavelength of incident light to mini- 
mize uncertainties in the scattering theory. The system is automated 
for data acquisition of the scattering matrix elements over an angular 
range from 2 degrees to 178 degrees of scattering. Both laser output and 
detector sensitivity are monitored to normalize the measured light scattering. 

An analysis of molecular scattering was performed to demonstrate the 
instrument's capabilities. The measured particulate angular scattering 
elements were generally found consistent with previous measurements. 

The primary need for improvement in the nephelometer design Is an 
increase in structural rigidity and a decrease in thermal expansion effects 
In order to achieve a more routine operational capability. However, the 
results of the calibration demonstrate that realistic measurements can be 
made with the nephelometer. Analysis of particulate measurements will 
appear in a subsequent paper. 

The authors would like to acknowledge Dale Hall, the electronics 
technician for the project. This research was funded by the Office of 
Naval Research under Grant Number N00014-75-C-0208. 
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I s = SCATTERED LIGHT 
lj = INCIDENT LIGHT 
0 = THE SCATTERING ANGLE 

Fig. 1. Scattering angles (where 0 is the azimuth 
angle of the scattered vector projected 
onto the transverse plane of the incident 
light and is measured from an arbitrarily 
chosen axis in the transverse plane). 
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Fig. 2. Total angular radiance scattering 
pattern for Mie particles for 
unpolarized incident light. 





Light scattering and detection schematic. 
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. Molecular scatter for parallel polarization of incident light. 
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Fig. 10. D 21 matrix element for the scattering 
volume. 
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